Mon. Not. R. Astron. Soc. 000, (0000) Printed 5 February 2008 (MN M^X style file v2.2) 



Constraining the Location of Microlensing Objects 
towards the LMC through Parallax Measurement in 
EAGLE Observations 



O ■ T. Sumi, 1 and Y. Kan-ya, 2 

■ 1 Princeton University Observatory, Princeton, NJ 08544-1001, USA; e-mail; suraiiSastro.princeton.edu 

' 2 Division of Earth Rotation, National Astronomical Observatory of Japan, Tokyo 181-8588, Japan; e-maihyukitoshi. kan-ya@nao.ac.jp 

bJO' 
< 



(N 



> 

ON 
ON 

m 
oo 
O 

(N 
O 

Oh 
6 



13 



Accepted Received in original form 



ABSTRACT 

We investigate the possibility of determining whether microlensing objects towards 
the Large Magellanic Cloud (LMC) are in a Galactic thick disc, or are in a Galactic 
halo, by using parallax measurements with an Earth-radius scale baseline. Our method 
makes use of EAGLE (Extremely Amplified Gravitational LEnsing) events which are 
microlensing events with an invisible faint source. We show that the rate of EAGLE 
events is as high as that of normal microlensing events, even if they are caused by 
dark stars in the Galactic thick disc. We explore the possibility of measuring the 
parallax effect in EAGLE events towards the LMC by using the Hubble Space Telescope 
(HST) or the Very Large Telescope (VLT). We find that EAGLE events enlarge the 
opportunity of parallax measurements by 4 ~ 10 times relative to that in normal 
microlensing events. We show that the parallax effect can be measured in ~ 75% 
(from the HST) and - 60% (from the VLT) of all EAGLE events if most lenses are 
stars in the Galactic thick or thin disc, while ~ 20% (from the HST) and ~ 10% (from 
the VLT) can be measured if most lenses are halo MACHOs. In combination with the 
finite source size effect observations, we can strongly constrain the location of lenses. 
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1 INTRODUCTION 

Several groups have carried out gravitational microlensing 
observations towards the Large Magellanic Cloud (LMC) in 
order to search for MAssive Compact Halo Objects (MA- 
CHOs) in the Galactic halo. Until now, 13-17 candidates 
have been found towards the LMC and the microl ensing op- 
tical dept h t from the events is 1.2±g;| x 1(T 7 ( |Alcock etj 
al. 20001 ) . The estimated typical lens mass depends on the 
adopted Galactic kinematic model ranging ov er 0.01 — 1 Mq 
( |Alcock et al. 2000b| ; [Honma fe Kan-ya 199^ ). 



Here M is the lens mass, x = Dd/D s is the normalized 
lens distance and Dd and D a are the observer-lens and the 
observer-source star distances. D s is hereafter assumed to 
be 50 kpc. The time variation of the parameter u = u(t) is 



u(t) 



t-t 



(3) 



We have learned that there are lens objects along the 
line of sight towards the LMC. However, the issues of where 
lens objects are and what they are, are still unclear. This is 
because a degeneracy occurs in ordinary mi crolensing event s 
for which the amplification is described by ( Paczyhski 198C ) 



where /3, to, £e = Re(M, a;)/i>t and v t are the minimum im- 
pact parameter in units of Re, the time of maximum mag- 
nification, the event time-scale and the transverse velocity 
of the lens relative to the line of sight towards the source 
star, respectively. From a light curve, one can determine the 
value of 0, to and t-E, where M, x and v t are degenerate in 
tE- This three- fold degeneracy is the essential difficulty in 
determining the nature of the lens objects. There are o nly 
marginally possible candidates, viz. ol d white dwarfs (Al- 



— for «< 1, 



ill 



cock et al. 1997, 2000b; Hansen 199£), old brown d warfs 



(lfio*rmia"'fc*l^an^va 1990 and primordial black holes ( [oka. 



where I, is the projected separation of the source and lens 



Tanaka fc Nakamura 200C ) 



in units of the the Einstein radius Re, which is given by 



Re{M,W) 



4(7 M n ... 



-xy. 



VT 



Possibilities for non-halo lensing objects have also been 
discussed, for example, d ark objects in a dark heavy com- 
ponent of the LMC itself (jAubourg et al. 1999| ; |Gyuk, Dalai 
fc Gricst 2000|; [Alcock et al. 200l|). There is also the possi- 



2 T. Sumi and Y. Kan-ya 



bility that the lenses are dark objects in the Galactic thick 
disc. While the microlensing optical depth r by known pop- 
ulations of sta rs in the Galactic th in disc and thick disc is 
of order 10 -9 (Alcock et al. 2000b), a maximal heavy thick 



disc, w hich may be surrounded by an extended dark halo 



comp os jed of particles, is also a possible Galactic compo- 
nent as the reservoir of lenses ( |Gatcs, Gyuk fc Turner 1995] 
Gates, et al. 1998). Such a thick disc can have r — 7 x 10~° 



( Pould 1994a| ; pould, Miralda-Escude fc Bahcall 1994| ), so 



the summed optical depth including the contribution from 
the Galactic thin disc (~ 2 x 10~ 8 ) and the LMC itself 



(~ 1 x 10 , 3ahu 1994 ) can be close to the observed value. 

The three-fold degeneracy can be resolved in som e kinds 
of exotic microlen sing events, e.g., the binary event ( Hardy 
Warme r 19951 lAlbrow et al. 1998|; |Alcock et al. 1999a| 



Honma. 1991% " Rionso et al 



20001; |An fc Gould 200 1|) and 
the finite source transit event (pould 1992|, |1994a|; pSfemiroflj 
fc Wicfrramasinghe 1994] [Witt fc Mao 1994 l p e n g 1997| ) . 
Sumi & Honma (200C ) pointed out that an extensive transit 



events search would make it possible to discriminate between 
the lenses in the Galactic halo and in the LMC. 

The third example is an event with parallax effect, 
which is essentially detected through the difference of light 
curves due to the spatial shift of an observer or observers. 
In such events we could determine the "reduced transverse 
velocity" v = vt/(l — x) of the lens. Examples of the paral- 
lax effects are the following. The annual parallax effect due 
to the Earth's motion around the Sun during an event leads 
an asy mmetry in the light curv e (e.g. pould 199 2j; pould 



Mirald4-Escude fc Bahcall 1994 Miyamoto fc Yoshii 1995] 



Alcock et al. 1995] |Mao 1999|; [Bennett et al. 20011; |Bond et 



al. 2001; iSoszyhski et al. 2001] ; pnhth, Mao fc Wozniak 2002] 
Mao et al. 2002). However, such events are rare. This paral- 
lax effect generally requires <e >~ 100 days, while <e ~ 40 
days for typical events. Another effect is parallax by the 
positional difference of two well-separated observation sites. 
In this case, we could measure the relative difference of the 
peak amplifications and the time at the peak amplifications 
between both. By observing an event from both a solar- 
orbit satellite and the Earth, we could utilize the parallax 
effect in almost every event, while it is difficult from two dis- 



tant observatories on the E arth (Refsdal 1966; Gould 1994b 



1995b; Holz fc Wald 1996). Furt 



rermore, the parallax could 
be measured by the positional shift of an observer due to 
the diurnal motion of the Earth, which was first advanced 
by An et al. (2002), and due to the orbital motion of an 



Earth-orbit sp ace telescope such as the Hubble Space Tele- 
scope (HST) (Honma 199!:) in binary events. In ordinary 
microlensing events this effect is quite small, but it is more 
efficient around the peak of high magnification events. 

For Galactic bulge events, Gould (1997) discussed the 
possibility of detecting the finite source size and parallax ef- 
fects by using two distinct ground-based telescopes in order 
to break the degeneracy in EME's (Extreme Microlensing 
Events, A > 200). With an approximate estimate for EME 
observations towards the LMC, he concluded th at it is not 
feasible . Howev er, with a more careful estimate, Nakamura 
& Nish] (1998|) showed that it is feasible to detect EAGLE 



(Extremely Amplified Gravitational LEnsing) towards the 
LMC. EAGLE is s imilar to the so-called "Pixel lensing" 
events (Gould 1996) but more simply defined as the events in 
which the source star is dimmer than observational limiting 



magnitude (ex. Vobs = 21 ~ 22), and not concerned whether 
the source star is resolved or not. Some EAGLEs would be 
EME's. EAGLE events could be efficiently detected with 
the "image s ubtraction method" or "Difference Image Anal- 



ysis (DIA ) " (|Alard fc Lupton 1998]; |Alard 
al. 19994 |20004 |Wozniak 200q ; 



2000 



Alcock et 



Bond et al. 2001), which 



has been recently developed and can perform more accurate 
photometry than DoPHOT and Pixel lensing method. So, 
we refer the term EAGLE in this paper. 

If lens objects are in the thin or thick disc (disc events) , 
Re projected onto the observer plane from the source star, 



Re(M,x) 



(4) 



is much smaller than that for halo events. In this case the 
light curve is more sensitive to the small displacement of 
the observer position. Therefore the fraction of parallax- 
measurable EAGLE events out of all EAGLE events for disc 
events will be much larger than that for halo events. This 
fraction is useful to discriminate statistically whether the 
lens objects are mainly in the thick disc or not. 

Here, we estimate the rate of parallax-measurable EA- 
GLE events towards the LMC. In §^ we summarize the ba- 
sic equations of EAGLE events. In EAGLE event rates 
are estimated. In §^ and ^ we describe the measurement 
of parallax effect from space and ground telescopes, respec- 
tively. In § ^ we calculate the fraction of parallax- measurable 
events. Discussion and conclusion are given in §0. 



2 BASIC FORMULAE FOR EAGLE EVENTS 

An EAGLE event is a microlensing event in which a source 
star is fainter than the observational limiting magnitude. A 
highl y amplified faint source w ould be detected as a new 
star. |Nakamura fc Nishi [1998 ) showed that the EAGLE 
event rate is fairly high towards the LMC. However an 
EAGLE search has not been involved in the microlensing 
surveys based on a DoPHOT-type PSF fitting photometry, 
which only measure already detected stars. 

A new CCD photometry method called DIA, in which 
an exposure frame is directly compared with a reference 
frame, enables much more accurate photometry at any place 
where any star isn't identified on the reference. This is more 
powerful for detecting EAGLE events than the DoPHOT. 
Many EAGLE events are expected to be found by DIA. 



2.1 Detection threshold 

A dim invisible source star with V-band apparent magni- 
tude V must be amplified brighter than EAGLE detection 
threshold Kh, which is slightly brighter than the observa- 
tional limit V^bs- So, the threshold amplitude is written as 



A T (Y) = 10' 



0.4(V-V th ) 



(5) 



A corresponding threshold impact parameter ut = 
ut{At(V)) which is the largest impact parameter to be de- 
tected as an EAGLE event de pends on V and V t ^. Thi s can 
be approximately written as ( |Nakamura fc Nishi 1998 ) 

-0.4(V-V th ) 



u T (V) ~ 10" 



(6) 
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3 GALACTIC MODEL AND EVENT RATE 

Here we give adopted models of the Galaxy and source stars, 
and estimate the event rate. Since only the relative event 
rate is discussed, common constant factors C, /, and the 
normalization of the luminosity function are not essential. 

When we consider microlensing by lens objects in the 
maximal disc ('disc events'), we adopt an exponential disc as 
the mass density distribution for the Galactic thin and thick 
disc. Since the line of sight to the LMC is almost tangential 
to the azimuthal direction in the Galaxy, the mass density at 
the region efficient to microlensing is not strongly dependent 
on the Galactocentric radius r but only on the height from 
the disc plane z, as long as the disc scale height is much 
smaller than r ~ 8.5 kpc. So we assume the disc density 
distribution only depends on z as 



§ exp H) 



(7) 



where E is the local disc column density and h is the 
scale height of the disc. For the thin and thick disc we 
set /ithin = 350 pc and h thi ck = 1400 pc respective ly. A 

2 (cf. [KuT 



certain acceptable va lue is E t hm — 50 M© pc 
jken fc Gilmore 1989| ). If we adopt this value, due to the 
maximal disc limit of ~ f00 M© pc~ 2 estimated from the 



rotation speed of the Galaxy (Binney fc Tremaine 1987 ) 
2 flGilmorc fc Reid 198$ pould 1994a 



50 M, 



Gould, Miralda-Escude fc Bahcall 1994| ). However, the value 



of t he local column density for each disc may be different 
(cf. iKui.jken fc Gilmore 1989!; |Bahcall, Flynn fc Gould 1992 



Creze et al. 1995 ; Holmberg & Flynn 2000). So we take E as 



a model parameter assuming that the combined mass of the 
thin and thick discs does not exceed the maximal disc limit. 
We change the fraction of each component within this limit, 
i.e., E thi n(E thick ) = 30(70), 50(50) and 70(30) M p C - 2 . 

For the disc events we assume the power-law mass func- 
tio n 4>(M) of lens objects defined by equations (9) and (10) 



Sumi & Honma (2000). The number density of lens ob- 



jects with the mass between M and M + dM is given by 
n(M,x)dM = f p{z(x))(j>(M)dM / M , where / is the mass 
fraction of the lens objects to the total mass and is assumed 
to be constant. We take the upper limit of mass M u = 50 
Mq and treat the lower limit M\ and the power-law index 
as model parameters as M\ = 0.1 or 0.01M© and ay = 2.35 
or 5, i.e., we assume both the ordinary (Salpeter IMF) and 
an extreme model for a darker population of stars as a con- 
stituent of the maximal disc. Our main results do not depend 
on this index very much (see § ^) . For the luminosity function 
of the source sta rs in the LMC, (jnJV) , we follow equations 
(15) and (16) in |5umi fc Honma (2000[ ) with the IMF index 
q s — 2.35, which is consistent with the l uminosity functi on 
of the LMC observed by using the HST ( [Holtzman 1997| ). 

For MACHOs ('halo events '), we adopt the sp herical 
'standard' halo model given by (Alcock et al. 2000b) 



2 1 2 
a +r 



(8) 



Phalo(r) = PO-5- 2> 

a 4 + r z 

where po = 0.0079 Mq pc -3 is the local mass density, r is 
the Galactocentric radius, ro = 8.5 kpc is the Galactocentric 
distance of the Sun, and a = 5 kpc is the core radius. We 
adopt delta-function mass functions with M = 0.1 and 0.5 
Mq because the lens mass function is not well known. 

For the observational parameters we assume Vobs = 21, 



which is a typical value for current microlensing programs 
with a 1-m class telescope. The events with the source star 
having V < 21 are normal microlensing events, while ones 
with V > 21 are EAGLE events. We consider two different 
values for the threshold magnitude; Vth = 19 and 20. 

The event rates of normal microlensing events (Fn) and 
EAGLEs (Fe), which are proportional to Re and utRe re- 
spectively, are given by 



T N = C / D s dx / dV 
Jo Jv, 



dMREV t MV)n(M,x), (9) 



r B = C/ D s dx 



dV 



dM«T RvVt </>l( V)n(M, a;), (10) 



where V u = 30 and V\ = 16 are the upper and lower limit 
of the luminosity function. Note that the integration is per- 
formed V < V \,s for Tn and V > Vobs for F3. 

We calculate the relative EAGLE event rate Ie/In for 
disc events, with the finite source size effect included. The 
results for Vth = 19 and 20 are r E /T N = 0.73 and 1.83, 
respectively, only weakly dependent on parameters of the 
disc structure (E,/i) and the mass function (ay, Mi). This 
ratio depends on the luminosity function of source stars. 
As a result, only in the most extreme case that Vth = 19, 
Qd = 5.0 and Mi = 0.01 M©, this ratio is slightly decreased 
due to finite source size effects. In the case that Vth = 20 
and ad — 2.35, we found that this effect is not negligible 
only when Mi < 10 -4 . For halo events, the events affected 
by finite source siz e effects are only seve ral percent out of 
all EAGLE events (Sumi & Honma 200C). So hereafter we 



neglect the finite source size effect. 

For events in which the source star is V > 25, the period 
in which the source is visible (V < 21) is less than 2 days. 
Then the detection efficiency for such events should be very 
low. So we estimated Fe /Fn in the case that the source star is 
V < 25, and found Fe/Fn = 0.59 and 1.47 with Vth = 19 and 
20, respectively. We checked this ratio in the conservative 
case of a gentler slope a s = 2.0 with Vth = 20 and V < 25, 
and found it to be 1.15. This is still sufficiently high. 

In this paper, we assume that microlensing events (nor- 
mal and EAGLE) would be detected by a 1-m class alert 
telescope and the real-time analyses with DIA, and then 
follow-up observations would be performed by an Earth- 
orbit space telescope or a large ground-based telescope. 



4 PARALLAX FROM SPACE TELESCOPE 

Here we describe the parallax effect observed with an Earth- 
orbit space telescope such as the HST. In the followi ng anal- 



ysis, w e use the coordinate system which is used in Honma 
(1999] ), i.e., the origin is set to be the centre of the Earth, 
and z-axis is set to be in the direction of the source star. 
The x-axis is set to be perpendicular to both of the z-axis 
and the orbital axis of the space telescope. The inclination 
i of the telescope orbit is defined as an angle between the 
z-axis and the orbital axis. We also assume that the space 
telescope is in a circular orbit with radius r st and angular 
velocity ui. The position of the telescope projected onto the 
observer plane (x-y plane) is written as 



T = 



[r s t cos(u;t + 8),r s t sin(urf + 5) 



cos i]. 



(11) 
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Here the angle 8 describes the position of the space 
telescope at t — to- Observed from the centre of the Earth, 
the position of the lens object on the lens plane is given by 

Le = [vt{t — to) cos 9 — b sin 9, v% (t — to) sin 9 + fecos 9], (12) 

where 9 is the angle between the x-axis and the direction of 
the transverse velocity of the lens vt and b is the minimum 
physical distance between the lens and the source-Earth line 
of sight. The position of the lens relative to the space tele- 
scope on the lens plane is given by 

L T = L E - (1 - as)T. (13) 

Then from the space telescope u(t) is written as 

f 1 ~l 2 

u'(t) — |~(* — *o) cos6> — /3sin# — £Cos(ut + 5)j 

+ |^-(* — to)sin# + {3cos9 — esin(wi + 5)cosi| , (14) 

where e = r st /7?E- Equation ( |l4| ) shows that the parallax 
effect due to the telescope motion causes a wavy trajectory 
of the lens object relative to that from the centre of the 
Earth. For orbital parameters for the space telescope, we 
assumed r st = 7000 km, i — 30° and P or b = 2h/oj — 97 min, 
which are similar to those of the HST. 

We show sample light curves without (A(«),thin line) 
and with (A(u'), thick line) the parallax in Fig. hp,, as well 
as the difference of the two 8 A = A(u) — A(u) in Fig. ^b 
in the case with tE = 5 days, vt = 30 kms -1 , x = 0.01, 
(3 = 0.05, 9 = 90° and 5 = 0°. Hereafter in the case using 
the HST, we take 5 = 0° since 8 is a random parameter 
related with just a phase of orbit. We ensured that this does 
not affect our main results. In the light curve from the space 
telescope, the small wavy perturbation from that observed 
from the centre of the Earth can be seen. In the SA light 
curve in Fig. Qb, we can clearly see that the fluctuation due 
to the parallax is enhanced around the peak. 

From the observation, we can derive the e and 9 in ad- 
dition to other parameters £e, /3, to and fo- From e and ts 
we can derive the "reduced velocity" v as follows 



(15) 



The reduced transverse velocity represents the projected rel- 
ative transverse velocity between the source star and the 
lensing object and will be of great use in investigating the 
lens location because this value for each component (halo, 
thick disc, and thin disc) is different (about 240 kms -1 , 
50 kms -1 and 30 kms -1 , respectively) and hence for disc 
events v is ~ 8 times smaller than that for halo events. 

The rate of parallax-measurable microlensing events de- 
pends on the uncertainty of v which is given by 



a 5 



/ dv\ 2 ( dv 

{'•all + r-aS 



(16) 



where <r e and <Jt E are the errors in the e and <e. We assume 
that the required threshold accuracy to measure v due to 
the parallax effect is 50%, i.e., cr^/v < 1/2, in this paper. 




t-t„ (day) 

Figure 1. (a) Light curves observed from the centre of the Earth 
(thin line) and from the HST (thick line), (b) 8 A light curve, (c) 
and (d) show integrands of the b tt and 6 E g, i.e. db ee /dt in 10 10 
e — day -1 and db t g/dt in 10 7 e _ day — 1 respectively, in the case of 
t E = 5 days, v t = 30 kms" 1 , x = 0.01, /3 = 0.05, 9 = 90° and 
8 = 0°. The amplitude of the modulation depends on e. 



The general discussion to estimate the error or the vari- 
ance of a set of para meters a, in fit t ing a distribution F(t, o,i) 
is presented by, e.g., Gould (1995a, 1998) using the minimum 
variance bound given from a well-known theorem in statis- 
tics. The covariance matrix with respect to for a series 
of measurements F(tk) at time with error cr/, is given by 



^ _ 2 dF(tk)dF(t k ) 



den 



do 



(17) 



The variance of at is just the diagonal elements cu. 

We take F(t) as the number of the detected photo- 
electrons at time t from the source star. During a short time 
interval T, we get F(t) = (fo(V)A(t; (3, t E , t ) + f b )T, where 
fo denotes the average photo-electron flux from the unam- 
plified source star and /{, denotes the background flux in the 
PSF aperture from the sky and the unlensed blending stars. 
A(t; P, tE, to) is the amplification produced by microlensing 
given by equations ([j]) and (|^). By using the DIA, the sig- 
nal is expressed as AF(t) = fo(A — 1) and (Tpit) = \J F(t). 
Taking the limit T — + 0, bij in equation (^) is given by 



fo(V)A + f b 



dfo(A 



dai 



a/ (A-i) 



dt. (18) 
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A parallax effect tends to be measured accurately in the case 
that /o is large, and/or j3, vt, x and/or M are small. 

We show the integrands of b tt and b t e in Fig. and [j]d, 
respectively, in the case of the source magnitude of V = 21 
and with the same parameters as in Fig. ^a. Here we as- 
sumed fo(V = 20) = 350 e~ s" 1 and ft, = 10 e~ s _1 , which 
is similar to those of the HST with the Advanced Camera 



for Surveys (ACS) (Pavlovsky 2001). ft, in the aperture of 
the PSF is estimated by taking account of the following 
sources; (1) Zodical light Vzl = 23.3 magarcsec -2 , which 
is the smallest because the Ecliptic latitude of the LMC is 
~ 90°. (2) Earth shine from the limb of the sunlit Earth 
Ves = 21.4 magarcsec -2 , which is the mean when the an- 
gle #n m b between the target and the bright Earth limb is 
larger than 12°. (3) Light from Blending stars Vlmc = 22.01 
magarcsec~ 2 , which is the mean V-band surface brigh tness 
of the inner 10 deg 2 of the LMC ( |de Vaucouleurs 1957| ). We 
reduced fo(V) and /;, by a factor of 3/4 because the LMC 
could not be observed when the HST is in shadow of the 
Earth and we limit 8nmb > 12°. (Note: the phase when the 
observation would be interrupted is not always the same, 
because the actual orbit of the HST is not circle but the 
latitude of the HST is changing between —30° and 30°.) We 
multiplied fo{V) and ft by 0.9 taking account of the dead 
time of camera assuming ~ 5 minutes exposures. 

In Fig. [J and [t|d, we can see that most of the informa- 
tion of e is in very short period (~ £e/3) around the peak, 
and the correlation between e and 8, i.e., b e e is negligible in 
a day-period observation. Furthermore, we ensured that e is 
completely independent of other parameters. 

We assume that the follow-up observation for EAGLE 
events would be carried out by the HST for one day from 
to- The extensive follow-up observations from ground-based 
small telescopes around the world make it possible to know 
to. We evaluate the uncertainty of the parameters a ai (a» = 
e, 8, £e). In Fig. ^, we show cr ai /oj as a function of /3 for the 
case of the typical event time-scale of = 40 days with 
same parameters used above except for 8 = 45°. The a f 
is the minimum at 9 = 90° and the maximum at 8 = 0°. 
Hereafter, we fixed 8 = 45° since 8 is a random parameter 
We ensured that the difference of the main results in § 6.z 
between the case using fixed 8 — 45° and the case taking 8 
at random, is negligible. 

Meanwhile £e is highly correlated with other parame- 
ters, i.e., f3, to and fo- If the event is observed symmetrically 
in time around to, £e is completely independent of to because 
to is an odd parameter in time. However, in the case that the 
event would be observed for a day from to, is no longer 
independent of to- In Fig. |2| we show (i) a tE /t-E (dotted line) 
and 0"t E /tE under the condition that (ii) fo (3dotted-dashed 
line) or (iii) both fo and to (dot-dashed line) are externally 
constrained. This (i) <Jt E in which there is a fo would be 
much improved over using two different telescopes to mea- 
sure the parallax in which there are two fo- Though, com- 
paring (i) and (ii), the former one deteriorates relative to 
the latter one because of degeneracies among the parame- 
ters t.E, (5 and fo, which is severe i n the case that the even t 
is observed only around the peak (Gould 1997, Han 1997). 




This can be improved by constraining fo by follo w-up obser- 
vation with the HST after the event (Han 1997). Hereafter 
we assume that fo would be constrained by the follow-up 
observation. As shown in Fig. |^, case (ii) is larger than (iii), 



Figure 2. Uncertainty cr ai /ai as a function of the impact param- 
eter P for e (solid line), 8 (dashed line), (i) tg (dotted line) and 
tE under the condition that (ii) fo (3dotted-dashed line) or (iii) 
both fo and to (dot-dashed line) are externally constrained, in 
the case of t B = 40 days, x = 0.01, Vt = 30 kms -1 , 8 = 45° and 
the source magnitude V = 21. The vertical dashed line indicates 
the EAGLE detection threshold impact parameter u^. 



but still quite small relative to cr e /e. Especially in the region 
of cr e /e < 1/2, where we are concerned, a tE /tE is at least 
10 times smaller than a c /e. In the case for brighter source 
event (V < 21), the contribution of at E /tE becomes slightly 
larger at a e /e ~ 1/2, but in this case to could be constrained 
well from the overall light curve taken by ground-based tele- 
scopes. Then hereafter we neglect a tE /tE in equation (|l6|), 
i.e., we can rewrite equation (hsh to a^/v ~ cr e /e. 

From the curve of cr e /e, we can obtain the critical im- 
pact parameter /3 cr i t to detect an parallax effect with 50% 
accuracy, which corresponds to a e /e = 1/2. 



5 PARALLAX FROM GROUND TELESCOPE 

Here we consider parallax measurements from an 8-m class 
telescope at a latitude of —30° on the Earth such as the Very 
Large Telescope ( VLT) . In this case the light curve would be 
expressed with the same equation as that from the HST (see 
equation (|l^)) except that the radius of the orbit is i? rb = 
R® cos 30° = 5500 km, the period is P or b = 2n/uj = 1 day 
and we can observe events only during night-time. Here R^ 
is the Earth radius. 

In Fig. |§| we show sample light curves of 5 A = A(u') — 
A(u) (a and d) and integrands of b te (b and e) and 6 e/ 3 (c and 
f) in the case of tE = 40 days, v t = 30 kms -1 , x — 0.01 and 
P = 0.05. Here we assumed that observations are performed 
for 7 hours a day. Fig. |l^.-c and [jji-f represent the case that 
the dark-side of the Earth is near to the lens relative to 
the centre of the Earth, i.e. #n — 9 — 90°, and the case 
that #n = 8 respectively, where #n is the angle between the 
centre of the observation night relative to the centre of the 
Earth and the x-axis. In this case, we can observe only some 
parts of the wavy light curve, and then e is highly correlated 
with other parameters, i.e., 8, tE, 13, to and fo- In the case of 
Fig. ^a-c, b e fs becomes large as well as 6 ee ,then <r s becomes 
significantly larger than the case of Fig. pp-f in which b^p 
becomes quite small, as well as the other cross terms bij. 
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Figure 3. Light curves of SA (a and d) in 10 -2 and integrands 
of b ee (b and e) and b € p (c and f) in 10 -11 in the case of tj; = 
40 days, v% = 30 kms -1 , x = 0.01, f3 = 0.05 and the source 
magnitude V = 21. Panels (a)-(c) and (d)-(f) represent the cases 
of 6jv — 9 = 90° , and the case of 0jv = & respectively. 



We estimate cr E by the same procedure in §|4|, i.e. eval- 
uating the covariant matrix. In this case, we assume the 
photo-electron flux is fo(V = 20) ~ 2000 e~ s" 1 using an 
8-m class telescope and thin CCD cameras with a filter near 
the V band with a relatively narrow bandwidth (~ 100 nm). 
The very broad-band filter is not advisable because that 
might not be able to eliminate the effect of different ial re- 



frac tion or differential extinction of the atmosphere ( Glould 
199 8j). | Ve reduce fo(V) and /(, by a factor of 2 to take ac- 
count of weather conditions and CCD camera dead time. We 
adopt the background flux /b = 1200 e~ s _1 corresponding 
to V = 20.6 mag, which is estimated as follows: the mean 
V-band surface brightness of the inner 10 deg 2 of the LMC 
is V = 22.01 magarcsec -2 , and the sky value is V = 21.6 
magarcsec -2 (ie Vaucouleurs 1957). We assume the aper- 
ture of the PSF is 0.49-7T arcsec (~ 0.7" seeing) which is 
similar to the typical seeing of the VLT. The total bright- 
ness in the aperture is 20.6 mag. We adopt this value for /b. 
Here we assumed that the follow-up observation would be 
taken for three nights ( 7 hours a night) around the peak. 

The LMC is visible enough (~ 7 hours) only during the 
southern summer (between equinoxes) from the site at a lat- 
itude of —30°, which corresponds to that #n is distributed at 
random between 0° ~ 180°. <r e is maximum at #n = 0° and 



180° and minimum at #n — 90° because of the inclination 
of the orbit axis i — 30° . Hereafter we use the mean value of 
#n — 45° . We ensure that the difference of the main results 
in the following analyses between the case using the fixed 
#n = 45° and that using a random #n is negligible. Anyway, 
this effect is much smaller than the dependence on 9 — (?n- 

In Fig. ^| we show <r e /e under the condition that fo is 
externally constrained (thick line), <r e /e under the condition 
that all other parameters are externally constrained (thin 
line), which is shown to compare with the case using the 
HST, and at E /tE under the condition that fo is externally 
constrained (dashed line) as a function of f3 in the case of 
i E = 40 days, x = 0.01, v t = 30 kms" 1 , 9 N = 45°, 8 = 45° 
and the source magnitude of V — 21. Upper and lower panels 
represent the case of 9 — —45° (i.e. #n — = 90°, which 
correspond to Fig. |l|a-c) and the case of 9 = 45° (i.e. #n = 9, 
which correspond to Fig. |l|d-f), respectively. In this figure, 
az/v heavily depends on #n — 9, it especially deteriorates in 
the case of #n — 9 = 90° relative to the thin line. In this case 
ct E /tE is also negligible. 

S = 45° is the most optimal case, i.e. one observes the 
peak at #n in the second observation night. In actual S de- 
pends on random parameters to and #n, and how early one 
can start the follow-up. Fig. |H| show a e /e as a function of S 
in the case using same parameters in Fig. ^ and (5 = 0.04, 
for the source magnitude of V = 20, 21, 22, 23 and 24 (from 
bottom to top). Here we take the mean value of a € /e over 
—45° < 9 < 135° because the 9 at which a £ /e is the max- 
imum or minimum, is related with S. In Fig. || we can see 
<r e /e does not depend on 8 very much, as long as observa- 
tions are spread roughly around to, i.e. it does not need to 
cover the peak exactly. When the follow-up observation is 
delayed more than 1 day (8 <~ —300), it becomes worse es- 
pecially for dimmer source events. This is worse for shorter 
time-scale events with much smaller /3 because most of the 
information of the parameters is contained in a very short 
period (~ t^P) around the peak. Furthermore in extreme 
cases (£e <C 10 days, /3 <C 0.01), generally <r E /e is small, 
however, it depends on whether the peak is covered exactly 
or not. This effect is large in dimmer source events (V > 25). 

From the curve of a f /e we can also obtain /3 cr it- We take 
the mean of /3 cr it over —45° < 9 < 135°. Uncertainty <r e also 
depends on the direction of v t . Note that this is negligible 
in the case using the HST. We take the mean of /? C rit on 
the cases of —Vt and +vt (i.e. same and opposi te d irection 
relative to the Earth rotation, respectively) in 



6.2 



6 FRACTION OF PARALLAX-MEASURABLE 
EVENTS 

Here we calculate the fraction of parallax-measurable events 
out of all events in the case of V t h = 20 and l£bs = 21. We 
adopt the same Galaxy density model and the source lu- 
minosity function as §H. To make the estimation realistic, 
we applied a photometric error 17% la rger than the p hoton 
noise, which is estimated for the DIA (Wozniak 2000). Fur- 
thermore, we specify the typical (transverse) velocity of the 
lens objects to estimate parallax-measurable event rate. 
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Figure 4. Uncertainty cr € /e under the condition that fo (solid 
line) and all other parameters (dotted line) are externally con- 
strained, and 0t E /iE under the condition that /o is externally 
constrained (dot-dashed line) as a function of j3 in the case of 
t E = 40 days, x = 0.01, v t = 30 kms" 1 , 6» N = 45°, <5 = 45° and 
the source magnitude V = 21. Upper and lower panels represent 
the case of 9 = 45° (i.e. 9^ = 9, which correspond to Fig. |l]a-c) 
and the case of 9 = 135° (i.e. (9jm — 9 = 90° , which correspond to 
Fig. [lji-f) respectively. The vertical dashed line indicates ur-. 



6.1 Model kinematics 

For disc events the typical transverse velocity of a lens ob- 
ject vt is estimated as follows. First we assume that the thin 
and thick discs have the same rotation velocity and the rel- 
ative velocity between the mean flow of the lens objects and 
the local standard of rest (LSR) is zero. In fact the effect of 
the drift velocity of the thick disc relative to the thin disc 
due to the large velocity dispersion of the thick disc is neg- 
ligible compared to other effects. v t is the composition of 
the tangential component of the mean velocity of the Earth 
relative to the LSR and the velocity dispersion of the 
disc er, which is assumed to be isotropic for both discs, as 

v 2 t {z) = 2a 2 (z)+v 2 B . (19) 

We take = 30 kms" 1 , which includes the orbital velocity 
of Earth around the Sun of 30 kms -1 and a velocity of 
the Sun relative to LSR of 20km s" 1 (Miyamoto & Zhu 



199c ). a 2 (z) and v 2 are assumed to be a function of z only 
(see the discussion on equation ((?])). The density profile of 
the disc is same as in § pi i.e., the exponential disc density 
distributions in equation (^) for thin (/ithin = 350 pc) and 
thick (ftthick = 1400 pc) discs with variable values of the 
local disc column density as E thin (E th i ck ) = 30(70), 50(50), 
70(30) Mqpc -2 . Solving the Poisson equation for the density 
profile composed of the thin and the thick discs cr 2 (z) is 
obtained from the z-comp onent of the Jeans equation (e.g. 
Binney fc Tremaine 1987 ) as 



ttG 



Ethi " exp {-htj) + EthiA exp 



2 Pthin(z) + Pthick(z) 

For halo events we adopt the same model as in § ^ and 
we assume vt = 190 kms -1 or 220 kms -1 . vt = 190 kms -1 
is more plausible because of the slight contribution from the 
disc gravity to the rotation curve, which is assumed to be 
flat with Viot = 220 kms . 



(20) 




-100 -200 -300 -400 -500 

6 (deg) 



Figure 5. Uncertainty cr € /t as a function of 5 in the case of us- 
ing the same parameters as Fig. Q and (3 = 0.04, for the source 
magnitude of V = 20, 21, 22, 23 and 24 from bottom to top. Here 
smaller <5 means that observations are delayed relative to the op- 
timal case that the peak is observed at the centre of the second 
night (<5 = 45°). 



6.2 Results 

We estimate the relative parallax-measurable event rate. 
The parallax-measurable event rate is written as 
pi pV n pM n 

T P =C D s dx dV dM 

JO JVi J Mi 

p cx «{V,M,x,vt{z{x)))RKVtfc{V)n{M,x), (21) 

where C is the constant common to the expression of Tn 
(equation (|^)) and Fe (equation (px|)). And Va = lobs for the 
normal events and V\ = lobs for EAGLEs. f3 CI i t is estimated 
by same procedure presented in § ^] and ^ for each parameter. 

In Fig. ^| we show the relative event rate distributions 
of parallax-measurable events dTp/dV from the HST (up- 
per panel) and from the VLT in the case of S — 45° (lower- 
panel), normalized by Fn (for V < V b B ) or Te (for V > V b B ) 
for thin disc, thick disc and halo events, as a function of 
source magnitude V and for various combinations of M\ 
and ad- We note that the absolute value of the vertical 
axis is not important because these distributions are rela- 
tive ones. These distributions are to be compared with the 
relative event rate distribution for all events (the bold dotted 
line). The left side relative to the vertical dashed line cor- 
responds to normal microlensing events and the right side 
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Figure 6. Relative differential event rate distribution of parallax- 
measurable events dTp/dV from the HST (upper panel) and from 
the VLT in the case of 8 = 45° (lower panel), normalized by In 
(for V < V b s ) or Ie (for V > V b s ), as a function of source 
magnitude V and for various combinations of M\ and aj. By 
the VLT, for a comparison, we also show that distributions in 
the case of 8 = —367.5° (i.e. the case that the observation start 
at to) f° r two cases of the thin and thick discs with Mi = 0.01 
and ay = 2.35 as same line styles (lower one). We also show the 
relative dT^/dV and dT^/dV by the bold dotted line. 



is for EAGLE events. For disc events the greater ay or the 
smaller Mi, the higher dTp/dV is. From Fig. |^ we can see 
that the event rate of parallax-measurable EAGLE events 
(V > 21) is much higher than that of normal microlens- 
ing events (V < 21). We can also see that for halo events 
the fraction of parallax-measurable events is still quite low 
even in EAGLE events. From the VLT (lower panel) , we also 
show that distributions in the case of 5 = —367.5° (i.e. the 
case that the observation starts at to) for two cases of the 
thin and thick discs with Mi = 0.01 and ay = 2.35 as same 
line styles (lower one). From this figure, we can see that 
follow up observations hv the VLT should he started ~ 1 



Table 1 . Ratio of the parallax-measurable EAGLE event rate to 
the parallax-measurable normal event rate for disc events. 

ay the thin disc the thick disc 

Mi = O.1M Mi = 0.01Af Q Mi = O.IMq Mi = O.O1M 

HST 



2.35 
5.0 



7.30 
6.71 



4.70 
4.05 



8.56 
8.07 



6.03 
5.36 



VLT 



2.35 
5.0 



8.44 
7.90 



5.78 
5.11 



9.50 
9.22 



7.24 
6.59 



Table 2. Ratio of parallax-measurable event in EAGLE events 
to that in normal events for halo events. 



u t =190km s" 1 «; t =220km s" 1 

M = 0.1M Q M = O.5M M = O.1M M = O.5M 

HST 



7.65 



8.29 



7.42 



8.01 



VLT 



4.58 



6.17 



4.20 



5.76 



events enlarge the opportunity of parallax measurements. 
The number of parallax-measurable EAGLE events is 5 ~ 9 
(from the HST) and 6 ~ 10 (from the VLT) times larger 
than that of normal events for disc events. For halo events 
these are ~ 8 (from the HST) and 4 ~ 6 (from the VLT) 
times larger than that of normal events. We hereafter refer 
Lp as the parallax-measurable EAGLE event rate because 
most of the parallax-measurable events are EAGLE events. 

In Table ^ we show the fractions of parallax-measurable 
EAGLE events out of all EAGLE events Fp /Ye and corre- 
sponding mean event time-scales (tv) for disc events in the 
case of V < 25 mag. For E th ick(Sthin) = 30(70), 50(50) 
and 70(30) M pc" 2 , the optical depths are 3.84(2.43), 
6.41(1.74) and 8.97(1.04) in 10" 8 , respectively, if / = 1, 
although the common factor / is not essential so long as a 
relative event rate is discussed, as noted in §|^. From Table 
H we see that (t E ) is strongly affected by Mi. Mi = O.O1M 
s eems to be consistent with the observed value tE ~ 40 days 
( |Alcock et al. 2000b| ) for the thick disc with a d = 2.35. 
Tp/Te only depends weakly on E thin | thick and ayf]. We can 
see that the parallax effect can be measured in ~ 75% (from 
the HST) and ~ 60% (from the VLT) of EAGLE events if 
the lenses are stars in the thick disc. 

The fractions for the thin disc are ~ 10% larger than 
that for the thick disc, though the contribution of these com- 
ponents is small because the optical depth i s about 25% rel- 



day bef pre to to detect the parallax efficiently. Hereafter, we 



ative to that of the thick disc (Gould 1994a; Gould, Miralda- 
Escude fe Bahcall 1994j). Mi might be between 0.1 and 0.01 



assume 8 — 45° for the case using the VLT. 

We evaluate the ratio of the parallax-measurable EA- 
GLE event rate to the parallax-measurable normal event 
rate r P (V > V obs )/r P (V < V obs ) in Table Fl for disc events 
with Ethin, thick = 50 Mq pc -2 and in Table H for halo events 
with the lens mass of M = 0.1 and 0.5 Mq in the case of 
V < 25 mag. From Table | and [| it is clear that EAGLE 



In our disc model, vt is larger for the larger Sthick case, and 
so lp /Te tends to be smaller. However for the thick disc from 
the HST it is inverse, because lp oc Dt, which is large for small 
x and this effect is slightly larger than the former effect in this 
case. Anyway, these effects are quite small relative to the others 
we are concerned. 
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Table 4. Fraction of parallax-measurable events in all EAGLE 
events and mean time-scales for halo events with mass functions, 
S(M): the delta-function with mass M, G(M): the Gaussian and 
G(logM): the log-normal Gaussian with mean mass M. 



Mass 


M 


ft = 


190km s" 1 


ft = 


220km s" 1 


Func. 


(M e ) 


<*E> 


r P /r E 


<*E> 


r P /r E 


HST 


<5(M) 


0.1 


22.5 


0.306 


19.4 


0.289 




0.5 


50.3 


0.201 


43.4 


0.191 


G(M) 


0.1 


46.0 


0.221 


39.8 


0.210 




0.5 


57.1 


0.194 


49.3 


0.184 


G{logM) 


0.1 


46.0 


0.251 


39.8 


0.238 




0.5 


51.2 


0.200 


44.2 


0.190 


VLT 


S(M) 


0.1 


22.5 


0.122 


19.4 


0.104 




0.5 


50.3 


0.105 


43.4 


0.095 


G(M) 


0.1 


46.0 


0.108 


39.8 


0.096 




0.5 


57.1 


0.100 


49.3 


0.091 


G(logM) 


0.1 


46.0 


0.108 


39.8 


0.095 




0.5 


51.2 


0.104 


44.2 


0.094 



Note: The event time-scales (4e) are given in day. 



from comparing the event durations with the observed value 
t-E ~ 40 days. Then, the parallax effect can also be measured 
in ~ 75% (from the HST) and ~ 60% (from the VLT) of EA- 
GLE events for the thin disc events. In Table ^[ we also show 
that fractions in the case of 5 = —367.5° and E t Mck,thm = 50 
Mq pc -2 from the VLT. In this case the fractions are ~ 30% 
smaller than that in the case of 5 = 45° . 

For halo events, we adopt the mass function as a delta 
function S(M), a Gaussian G(M) and a log- normal Gaus- 
sian G(logM) distributions because the lens mass function 
is not well known. For S(M) we adopt M = 0.1 and O.5M0. 
For G(M) and G(logM) we take the mean of the mass 
M = 0.1 and 0.5 Mq. The variance is 0.4M Q for the G(M) 
and log(O.4M /M) for G{logM). Of course the mass func- 
tion is zero for negative M in G(M). The typical transverse 
velocity is v t = 190 or 220 kms' 1 . The estimated fractions 
in the case of V < 25 mag are shown in Table ^. The optical 
depth is 4.8 x 10 -7 if / = 1. From Table 0, we can see that 
the fraction of parallax-measurable events out of all EAGLE 
events is ~ 20% (from the HST) and ~ 10% (from the VLT). 

The fraction of parallax-measurable events depends on 
the lens location. Then we can also statistically discriminate 
whether the lenses are in the thick disc or halo, by using the 
parallax-measurable EAGLE event rate. 



7 DISCUSSION AND CONCLUSION 

We have seen that the EAGLE event rate is as high as that 
for normal events even for disc events towards the LMC. 
Since the period in which sources are visible (V < 21) is 
usually short (1 day ~ 40 days), the detection efficiency 
heavily depends on the observational frequency. The obser- 
vational programs currently undertaken by most groups are 



not adequate. Hourly monitoring with a 1-m class dedicated 
telescope and the real-time detection of EAGLE events by 
the DIA are required to issue alerts with a high detection 
efficiency. However, for the events in which the source star 
is V > 25 this period is less than 2 days, so it is difficult 
to detect. Thus we estimated Fe/Fn and Fp/Fe in the case 
that the source star is V < 25. These are decreased but still 
sufficiently high. Of course, the larger alert telescopes make 
it easier and faster to issue the alerts. 

Estimating the parallax-measurable event rate, we ad- 
vocate follow-up observations with a space telescope such as 
the HST and an 8-m class ground-based telescope such as 
the VLT. We have found that EAGLE events enlarge the 
opportunity of parallax measurements by 5 ~ 9 (from the 
HST) and 6 ~ 10 (from the VLT) times for disc events, and 
by ~ 8 (from the HST) and 4 ~ 6 (from the VLT) times for 
halo events relative to that in normal microlensing events. 
We have also found we can measure the parallax effect in 
~ 75% (from the HST) and ~ 60% (from the VLT) of EA- 
GLE events if the lenses are in the thick or thin disc, and 
in ~ 20% (from the HST) and ~ 10% (from the VLT) if the 
lenses are in the halo. Since v for halo objects are 5 ~ 8 times 
larger than that for disc stars, we can determine whether the 
lens objects are in the halo or discs for each event. We can 
also statistically constrain the lens locations by using the 
parallax-measurable EAGLE event rate. 

In follow-up observations from the HST, in this paper 
we assumed that the observations start just after the peak as 
the most conservative case. Of course, an observation around 
the peak is better than that just after the peak. However 
predicting to is not so easy for very faint source events. Ex- 
tensive follow-up observations from small ground-based tele- 
scopes around the world are needed to predict to and inform 
to the HST. Furthermore a flexible operating program of the 
HST for the alerts are required. If the alert telescope is at 
a latitude —30°, the alert can be issued for half of the year 
(Southern Summer), while at a latitude —44° (such as New 
Zealand) it can be done all year round. The total operation 
time of the HST would be several days per year. 

From the VLT, we assumed 8 = 45°. However for the 
events with faint source (V > 24) the time until to is short, 
and the beginning of the follow-up would tend to be delayed. 
A delay of ~ 1 day reduces the possibility of measuring the 
parallax by ~ 30% as shown in §^^. This observation can 
be done for only half of the year (Southern Summer). 

The true source flux fo is needed to measure the precise 
value of v in the light curve fitting. Then follow-up obser- 
vations by a high resolution telescope such as the HST are 
needed to get an accurate fo after the event. 

In short, a practical observation strategy would be to 
observe hourly with a 1-m class telescope and perform real- 
time analysis with DIA to issue alerts to world observatories 
and the HST or the VLT for follow-up observations. Then 
after the events fo should be measured by the HST. 

To demonstrate this specifically we estimate the num- 
ber of expected parallax-measurable events for the two cases 
that these are mainly halo events or disc events. In both 
cases, the thin disc events are included. For the typical pa- 
rameters « d = 2.35, Ethin = Ethick = 50Af Q pc~ 2 , Vth = 20, 
a detection efficiency of 50% and source stars of V < 25, one 
can expect to find ~13 EAGLE events from 3-year observa- 
tions of 11 square degrees of the LMC central region (as the 
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Table 3. Fraction of parallax-measurable events in all EAGLE events and mean time-scales for disc events. 



ay Sthick the thin disc the thick disc 





(M pc- 2 ) 


M 1 = 


0.1M© 


M\ = 


O.O1M 


Mi = 


O.1M 


Mi = 


O.O1M 






(tE) 


r P /r E 


(tE) 


r P /r E 


(tE) 


r P /r B 


(tE) 


r P /r E 


HST 


2.35 


30 


61.7 


0.659 


23.2 


0.867 


128.9 


0.513 


48.6 


0.745 
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0.656 


21.0 


U.oDO 


1 1 fi o 


U.010 


A A Q 


n 7/1 7 

U. 1 4 1 


2.35 


70 


49.7 


0.654 
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0.749 


5.0 


30 


41.3 
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13.0 


0.938 


86.3 
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5.0 


50 


37.2 
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11.7 


0.936 


79.5 
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25.0 


0.834 


5.0 


70 
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23.1 
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VLT (<5 


= 45°) 










2.35 


30 


61.7 


0.486 


23.2 


0.704 


128.9 


0.355 


48.6 


0.580 


2.35 


50 


55.6 


0.479 


21.0 


0.687 


118.8 


0.356 


44.8 


0.575 


2.35 


70 


49.7 


0.469 


18.7 


0.664 


109.6 


0.355 


41.3 


0.566 


5.0 


30 


41.3 


0.555 


13.0 


0.785 


86.3 


0.421 


27.2 


0.668 


5.0 


50 


37.2 


0.546 


11.7 


0.765 


79.5 


0.420 


25.0 


0.660 


5.0 


70 


33.3 


0.532 


10.5 


0.736 


73.4 


0.418 


23.1 


0.648 








VLT (8 = 


-367.5°) 










2.35 


50 


55.6 


0.359 


21.0 


0.451 


118.8 


0.291 


44.8 


0.420 


5.0 


50 


37.2 


0.393 


11.7 


0.484 


79.5 


0.336 


25.0 


0.465 



Note: The event time-scales (is) are given in day. 



MACHO collaboration does). In the case of follow-up from 
the VLT, the expected number would be half. In these 13 
EAGLE events, ~ 2 events (15%) are due to the stars in 
the thin disc and a further ~ 11 events are due to MACHOs 
or the dark stars in the thick disc. In considering these 11 
events, reasonable parameters are M = 0.1 or 0.5Mq except 
M = Q.IMq in the case of the 8- function for halo events, and 
Mi = O.OIMp) for disc even ts, to be consistent with (£e) ~40 
days ( Alcock et al. 2000b ). In this case, we will be able to 
measure v in ~ 10 (from the HST) or ~ 4 (from the VLT) 
events for disc events, ~ 4 (from the HST) or ~ 1 (from the 
VLT) event for halo events, which include thin disc events. 
We can constrain lens locations strongly based on the 3-year 
statistics of these observations. 

In conclusion, one could statistically discriminate 
whether the typical lens locations are in a thick disc or not, 
using parallax measurements even with ~ R® scale baseline. 
One could also distinguish whether the lenses are MACHOs 
or stars in the LMC itself through finite source size effect s 
measurements in EAGLE events ( Sumi fc Honma 2000 ). 
Therefore one could identify lens objects as halo MACHOs, 
dark stars in the Galactic thick disc, or stars in the LMC 
through these observations. 

A real-time alert system with DIA, has been introduced 
by the MOA collaboration^ from 2000 and by the OGLE 
collaboration^] from 2002. 



see )t.t.p://Bw.phvs.r.aTit.aT-biiTv.ar..Rz/~ph.vsih/a1e-rt./a1ert. .html 



xttp : //www. astrouw. edu.pl/~ogle/ogle3/ews/ews .html 
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